Thermal orientation of electron spins 
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It is shown that the spin orientation of free electrons occurs in low-symmetry semiconductor 
structures if only the electron gas is driven out of thermal equilibrium with the crystal lattice. The 
proposed mechanism of such a thermal orientation of electron spins is based on spin dependence of 
the electron-phonon interaction which tends to restore equilibrium. The microscopic theory of the 
effect is developed here for asymmetric (llO)-grown quantum wells where the electron gas heating 
leads to the spin orientation along the [110] axis in the quantum well plane. 

PACS numbers: 72.25.Rb, 73.63.Hs, 75.75.+a 



I. INTRODUCTION 



II. MICROSCOPIC MODEL 
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The spin-related phenomena in semiconductors have 
been attracting considerable attention since the discov- 
ery of the optical orientation of electron and nuclear spins 
and the basic mechanisms of spin relaxation Much 
effort in this field is currently focused on the development 
of novel methods of spin orientation of free carriers rang- 
ing from the spin Hall effect Q and spin-dependent tun- 
neling j| to the optical orientation by linearly polarized 
light [4|. Here we show that the spin orientation of free 
electrons can be achieved in semiconductor nanostruc- 
tures by simple electron gas heating. We demonstrate 
that, in nanostructures of sufficiently low space symme- 
try, the electron-phonon interaction tending to restore 
equilibrium between the electrons and the crystal lat- 
tice is spin-dependent and leads to the spin orientation 
of carriers. Such a thermal orientation of electron spins 
is considered here for asymmetric (HO)-grown quantum 
wells (QWs) where the electron gas heating leads to the 
spin orientation along the [110] axis in the QW plane. 

The possibility to achieve the spin polarization of car- 
riers S caused by the disturbance of thermal equilib- 
rium in asymmetric (llO)-grown QWs follows from sym- 
metry analysis. Indeed, such structures belong to the 
point group C s that contains only two symmetry ele- 
ments, namely, identity and a mirror plane perpendicu- 
lar to the x axis. Here we use the following coordinate 
frame: x||[110] and j/||[00l] are the in-plane axes, and 
« || [110] is the growth direction. Reflection by the mirror 
plane changes the sign of the y and z components of the 
spin axial vector S but does not modify the x compo- 
nent. Therefore, the spin component S x is an invariant 
in asymmetric (llO)-grown QWs suggesting that the spin 
polarization along the x axis can emerge if the electron 
gas is driven out of equilibrium by any means. Partic- 
ularly, in the case of disturbance of thermal equilibrium 
between the electrons and the crystal lattice, the spin 
orientation can phcnomcnologically be described by 
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The microscopic model of the thermal orientation of 
electron spins is caused by the energy relaxation of hot 
carriers and includes two stages which are illustrated in 
Fig. 1(a) and Fig. 1(b). 

In the first stage [Fig. 1(a)], the carriers lose a part 
of their kinetic energy by emitting phonons. Such en- 
ergy relaxation processes shown by curved arrows are 
spin-dependent [5|. In QW structures without an inver- 
sion center, spin-orbit interaction adds an asymmetric 
term to the probability of electron scattering by phonons 
which is linear in the wave vector components [6(. As is 
shown in the next Section, the dominant spin-dependent 
contribution to the probability of electron scattering in 
(llO)-grown QWs is proportional to a z (k x + k' x ), where 
a z is the Pauli matrix, k x and k' x are components of the 
initial and scattered electron wave vectors. Due to the 
spin-dependent asymmetry of the electron-phonon inter- 
action, electrons with the spin "up" (along the z axis) 
predominantly vacate the excited states with positive k x 
while electrons with the spin "down" vacate the excited 
states with negative k x . This leads to a nonequilibrium 
distribution where the spin-up hot carriers occupy mainly 
the left-hand branch of the dispersion curve (carriers with 
the opposite spin orientation have gone to the subband 
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where AT = T e — To, T e and To are the electron and 
lattice temperatures, respectively. 



FIG. 1: Microscopic origin of the electron gas spin orienta- 
tion due to the energy relaxation of hot carriers, (a) Spin- 
dependent asymmetry of the electron-phonon interaction fol- 
lowed by (b) spin precession in the effective magnetic field 
leads to the spin orientation of electrons along the x axis. 
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bottom) while the spin-down carriers occupy mainly the 
right-hand branch. 

In the second stage [Fig. 1(b)], a net spin orientation 
of the electron gas appears as a result of the subsequent 
spin precession of noncquilibrium carriers in the effective 
magnetic field induced by the spin-orbit interaction Q. 
The field has nonzero in-plane components in asymmet- 
rically growth QWs, e.g., due to the Rashba effect Q. 
Therefore, the spins of noncquilibrium carriers, directed 
along or opposite to the z axis after the phonon emission, 
will process in the effective field as shown in Fig. 1(b). 
Note that electrons with the initial spin —1/2 and wave 
vector k x > are acted upon by the effective field with 
the Larmor frequency 17k, while particles with the ini- 
tial spin +1/2 and the negative wave vector feel the field 
with the Larmor frequency il-k- The effective magnetic 
field induced by spin-orbit interaction is an odd function 
of the wave vector, therefore, f2-k = — ^k, arid the spins 
of particles with positive and negative values of k x rotate 
in opposite directions. This gives rise to the spin compo- 
nent S x > for all hot electrons in the subband and the 
net spin polarization of the electron gas. 

We note that the electron spin caused by the energy 
relaxation processes is nonequilibrium and, therefore, can 
be rotated by an external static magnetic field, similarly 
to the Hanle effect in optics. 



III. THEORY 

The spin-dependent asymmetry of the electron-phonon 
interaction can be obtained if one takes into account 
both k • p admixture of the valence-band states to the 
conduction-band wave function and the phonon-induccd 
intcrband coupling. To first order in the k ■ p theory, 
where k = (k x , k y ) is the in-plane wave vector, the elec- 
tron wave function in a (llO)-grown QW has the form 
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Here S and X', Y', Z' are the Bloch functions of the con- 
duction and valence bands at the T point of the Brillouin 
zone, ijj and v = (v x ,v y ,v z ) are the envelope spinors. 
The Bloch functions X', Y', and Z' in Eq. © are re- 
ferred to the cubic axes z'||[100], y'||[010], and z'||[001], 
respectively. The envelope spinors of the valence-band 
states are related to the conduction-band spinor by 
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where P cv = (S\p z i\Z'} is the interband matrix element 
of the momentum operator, mo is the free electron mass, 
Eg is the band gap, A so is the valence-band spin-orbit 
splitting, and a is the vector of the Pauli matrices. 

We consider the electron scattering by acoustic 
phonons due to the deformation-potential mechanism. In 



cubic noncentrosymmctric crystals such as zinc-blcndc- 
type semiconductors, the strain induces a coupling be- 
tween the conduction-band and valence-band states [3, 
9] . The matrix elements of such a coupling have the form 

Vs,X' = ^cvUy'z', Vs,Y' = ^cvUx'z', Vs,Z< = ^. C vU x ' y ', 

where S c „ is the interband constant of the deformation 
potential, and u a p are the strain tensor components used 
here in the primed coordinate system. It is the strain- 
induced interband coupling together with the spin-orbit 
splitting of the valence band that leads to spin-dependent 
asymmetry of the electron-phonon interaction. Taking 
into account k • p mixing given by Eqs. © and ([3]) and 
allowing for the interband coupling, we derive for the 
Hamiltonian of the electron-phonon interaction in (110)- 
grown quantum wells 



Ki-phon^k) = ^ c 2_^u aa + £E cv {(k x +k' x )x (4) 

a 

[o z (uzz - u xx )/2 - o-yUyz] + {k y + k' y )[a x u yz + <7 z u xy }} , 

where 5 C is the intraband constant of the deformation po- 
tential responsible for the dominant (spin-independent) 
part of the electron-phonon interaction, and £ is the co- 
efficient given by 
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In deriving Eq. ([4J, we have expressed the strain tensor 
components in the primed axes via those in the QW coor- 
dinate system xyz by using the equalities u x i y i = (u zz — 
u xx )/2, u x i z , = -{u xy +u yz )/\/2, u y i z < = [u xy -u yz ) I 

The probability of electron scattering is determined by 
squared matrix elements of the electron-phonon interac- 
tion (J?)). The dominant contribution to spin-dependent 
asymmetry of the electron scattering in quantum wells 
is given by terms proportional to the u zz component of 
the phonon-induced strain tensor. This is because the 
strain tensor components depend on the wave vector and 
polarization of the phonon involved, and the in-plane 
component of the phonon wave vector q\\ = |k — k'| is 
typically much smaller than the out-of-plane component 
1z ~ tt/o, where a is the QW width. Thus, the principle 
contribution to the scattering asymmetry in (llO)-grown 
QW structures is proportional to <r z (k x + k' x ). This term 
is taken into account in calculations which follow. Wc 
note that k-linear terms in the matrix elements of the 
electron-phonon interaction can also be obtained in the 
second order in the k • p theory, as was done in Ref. 
for the electron scattering by charge impurities. However, 
these terms lead to no essential contribution to asymme- 
try of the electron scattering by phonons nor the spin 
orientation caused by the energy relaxation processes. 

As is shown in Section [ill spin-dependent asymmetry 
of the scattering processes followed by the precession of 
electron spins in the effective magnetic field leads to spin 
orientation of the electron gas. We assume that the spin 
relaxation time of carriers is much longer than the ther- 
malization time controlled by electron-electron collisions 
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which is in turn much longer than the momentum relax- 
ation time t p , and Q^t p <C 1. Then, the spin generation 
rate is given by [H 



s = ^ t p [rj k x g k ] 



(6) 



where gk = Tr[erG(k)]/2, G(k) is the spin matrix de- 
scribing the carrier redistribution in k-space due to the 
scattering by phonons. Components of the matrix G(k) 
have the form (see, e.g., Ref. [§]) 
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sl.kl q,± 



{^Jkl^lkl /(£kl)[l - /(£k)] <5( £k - e kl ± hu>J 
- KiktkKitU'k /(£k)[l - /(£ki)] 6(s kl - £k ± ^ q )} , 

where Veliki is the matrix element of the electron scat- 
tering assisted by emission (+) or absorption (— ) of a 
phonon, /(ek) is the distribution function of carriers, 
t k = h 2 k 2 /(2m*) is the kinetic energy, and w q is the 
phonon frequency. Note that spin-orbit splitting of the 
energy spectrum is neglected in Eq. (JT]). 

In (HO)-grown QW structures, components of the Lar- 
mor frequency of the effective magnetic field have the 
form 



n u - h h xy k y , lyx k x , lzx k x ) 
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The parameter ^ zx is caused here by the lack of an in- 
version center in the host crystal, while "f xy and 7 yx are 
non-zero due to the QW asymmetry only. We assume 
that electrons obey the Boltzmann statistics and the lat- 
tice temperature is not very low, fc^To 3> hujq. Then, 
one derives 
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where p is the crystal density, N e = 2^ k /(e k ) is the 
two-dimensional electron concentration, and f(z) is the 
function of size quantization. 

In the steady state regime, when the electron tempera- 
ture T e and the lattice temperature To are held constant, 
the spin density S x is determined by balance between 
the spin generation and relaxation processes; S x = S x T Xl 
where T x is the spin relaxation time. In (OOl)-grown 
QWs, the time T x is given by [HI 

T; "-f7B^ ( «* ) + ,flL,li " (10) 

where the angle brackets mean averaging over directions 
of the wave vector. For the Boltzmann distribution, 
Eq. (fT0|) assumes the form 
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Finally, we obtain for the steady spin density 



S x = - 
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The estimation for the average electron spin gives 
S x /N e ~ 10~ 5 for the electron temperature T e = 100 K, 
the ratio AT/T e m 1, the quantum well width a = 100 A, 
and band parameters m* = 0.07 mo, ^ = 0.4 A, y/h = 
10 5 cm/s, S c = —8 cV, E cv = 3 eV [l( corresponding to 
GaAs-based QW structures. Thus, for the carrier den- 
sity N e = 10 11 cm -2 , the spin density S x is of the order 
of 10 6 cm -2 which is well above the experimental reso- 
lution. We also note that the modest estimated value of 
the spin polarization is due to the fact that the energy 
relaxation by acoustic phonons is ineffective. The spin 
polarization can considerably increase if optical phonons 
are involved in the energy relaxation processes. 
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